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We propose an interferometric method for statistically discriminating between nonorthogonal 
states in high dimensional Hilbert spaces for use in quantum information processing. The method 
is illustrated for the case of photon orbital angular momentum (0AM) states. These states belong 
to pairs of bases that are mutually unbiased on a sequence of two-dimensional subspaces of the full 
Hilbert space, but the vectors within the same basis are not necessarily orthogonal to each other. 
Over multiple trials, this method allows distinguishing 0AM eigenstates from superpositions of mul¬ 
tiple such eigenstates. Variations of the same method are then shown to be capable of preparing 
and detecting arbitrary linear combinations of states in Hilbert space. One further variation allows 
the construction of chains of states obeying recurrence relations on the Hilbert space itself, opening 
a new range of possibilities for more abstract information-coding algorithms to be carried out ex¬ 
perimentally in an simple manner. Among other applications, we show that this approach provides 
a simplified means of switching between pairs of high-dimensional mutually unbiased 0AM bases. 
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I. INTRODUCTION 

Quantum key distribution (QKD) aims to enable two 
agents, Alice and Bob, to generate a shared crypto¬ 
graphic key while preventing an unauthorized eavesdrop¬ 
per, Eve, from gaining significant information about the 
key without being revealed. In the context of optics, the 
most common variable used to encode the key is the pho¬ 
ton’s polarization. Polarization spans a two-dimensional 
space, so that it normally can encode only one bit of key 
per photon. In order to generate multiple bits per pho¬ 
ton, other variables may be enmloyed which span higher 
dimensional Hilbert spaces [i|, 0 . The most promising of 
these is the photon’s orbital angular momentum (0AM). 
The 0AM along the propagation axis is Lz = Ih, where 
the topological charge I is restricted to integer values, 
and in principle can be arbitrarily large. If a stream of 
photons is produced whose I values are allowed to range 
between —L and +1/, then each photon can encode up to 
log 2 ( 2 I/ + 1) bits of information in its 0AM. 

In 0, an approach was introduced for using 0AM to 
generate a secret encr^mtion key based on a novel en¬ 
tangled light source [j-[8| that produces output with an 
0AM spectrum whose absolute values are restricted to 
the Fibonacci sequence. Recall that the Fibonacci se¬ 
quence starts with initial values Fi = 1 and F 2 = 2, 
generating the rest of the sequence via the recurrence 
relation Fn +2 = ^n+i + ^n- These states are used to 
pump a nonlinear crystal, leading to spontaneous para¬ 
metric down conversion (SPDC), in which a small pro¬ 
portion of the input photons are split into two entangled 
output photons, called the signal and idler. After the 
crystal, any photons with non-Fibonacci values of 0AM 
are filtered out. The result is that angular momentum 


conservation and the Fibonacci recurrence relation force 
the two output photons to have 0AM values that are ad¬ 
jacent Fibonacci numbers (for example Fm and F^+i). 
These photons can then be used to generate a secret key; 
see d, [sl for more details on the procedure. 

The Fibonacci key distribution protocol requires the 
ability to distinguish between single 0AM eigenstates 
and superpositions of pairs of eigenstates. These two 
types of states belong to two mutually unbiased bases on 
subspaces of the full Hilbert space. The state discrimi¬ 
nation cannot be done unambiguously on a single trial, 
just as it can not be determined in a single trial whether 
a photon’s polarization is vertical or diagonal; however, 
over many trials a statistical picture of the outcomes can 
be built up in such a way that eavesdropping alters the 
outcome probability distributions in a detectable man¬ 
ner. In this manner, an eavesdropper can be revealed 
over multiple trials even though it may not be possible 
to identify errors on any individual trial. The details of 
how the eavesdropper revelation works may be found in 
0 . 

The Fibonacci protocol is an example of a more ab¬ 
stract approach to QKD in which, rather than sim¬ 
ply modulating between two spatial measurement bases, 
more abstract modulations are carried out, switching be¬ 
tween more general sets of states in Hilbert space. Unlike 
polarization-based protocols, the relevant unbiased bases 
in Hilbert space do not correspond directly to projections 
onto bases in physical space, allowing more freedom in 
the choice of bases used. By opening up a broader range 
of states to manipulate, this approach has promise to in¬ 
crease the key capacity per photon, to allow new meth¬ 
ods for safeguarding the security of the key, and in some 
cases to simplify experimental implementations. In order 
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to carry out this program, it is necessary to have a sim¬ 
ple means of producing and detecting general (possibly 
nonorthogonal) linear combinations of basis states. 

In this paper, we examine the issue of discriminating 
between non-orthogonal states formed from superposi¬ 
tions of 0AM eigenstates. We propose an interferometric 
method which allows the user to distinguish between two 
different superpositions or to distinguish superpositions 
from individual eigenstates. There is some probability of 
error in the state identification, due to the nonorthogo¬ 
nal nature of the states involved, but for the purposes of 
QKD this can be dealt with in the reconciliation stage, as 
detailed in [9| . More general discussions of the problem of 
non-orthogonal state discrimination with minimal error 
may be found in UMll. Another approach to the study 
of superposition states in 0AM space via optical filtering 
with spatial light modulators can be found in Bill. In 
approach was taken to high-dimensional time- 
bin-encoded states that is similar in spirit to what is 
done here, allowing measurements to be done in arbi¬ 
trary bases. 

We first display a solution to the problem of sorting 
pairwise superpositions of states. We then generalize the 
apparatus to allow the detection of arbitrary linear super¬ 
positions of 0AM states. It is further shown how slight 
changes could turn the same structure into a means of 
synthesizing arbitrary superposition states. This is sim¬ 
ilar to the case of conventional phase and polarization 
coding in telecom systems, where Mach-Zehnder inter¬ 
ferometers are used for both coding and decoding of the 
signal. In the current case, the type of tree-like structure 
used to generate and detect these superposition states 
will be referred to as a superposition generation and de- 
teetion tree (SGDT). The use of different SGDT’s allows 
the generalization of the Fibonacci protocol to similar 
protocols based on other linear recurrence relations. Fi¬ 
nally, we move from recurrence relations on the real num¬ 
bers to recurrence relations on Hilbert space, showing 
how to construct chains of states (as opposed to eigen¬ 
values) that obey arbitrary linear recurrence relations; 
this is done via a set of nested structures we will refer to 
as a recursive state generator (RSG). 

This paper is organized as follows. In section [IT] we re¬ 
view the Fibonacci protocol as a motivating example for 
what follows. In section uni we then discuss the prob¬ 
lem of nonorthogonal state discrimination and introduce 
SGDT’s for the states relevant to the Fibonacci proto¬ 
col, with generalizations that cover other types of super¬ 
position states in section nvi We then introduce RSG’s 
in section m for the preparation of superposition states, 
and discuss the idea of recurrent sets of states in Hilbert 
space. Finally, we briefly discuss conclusions in section 
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FIG. 1: (Color online) Schematic setup for generating cryp¬ 
tographic key with Fibonacci-valued 0AM states. The L and 
D detection stages respectively test for different 0AM eigen¬ 
states or different superposition states. Interference between 
scatterings in aperiodic Vogel spiral lead to an 0AM spectrum 
in the pump that is restricted to the Fibonacci sequence. The 
filters after the crystal remove any non-Fibonacci values that 
appear in the down conversion process. See for more de¬ 
tailed description of the apparatus. 


II. FIBONACCI KEY DISTRIBUTION 


To provide a motivating example of the need for dis¬ 
tinguishing between nonorthogonal superposition states 
within unbiased bases, we review the essential ideas of the 
Fibonacci protocol Q, which allows multiple key bits to 
be generated per photon. In this protocol, distinguishing 
such states is necessary both to determine the key and 
in order to provide security against eavesdropping. Se¬ 
curity of the protocol stems from the random switching 
between measurements in two bases (the D and L bases 
introduced below) that are mutually unbiased on each of 
a chain of two dimensional subspaces; but in one of these 
bases the two basis vectors are mutually nonorthogonal. 

The basic setup of Q is shown schematically in Fig. [TJ 
The source on the left consists of an aperiodic nano-array 
of scatterers in the form of a Vogel spiral , followed 
by a nonlinear down conversion crystal. The Vogel spiral 
causes the scattered waves from the different scattering 
centers to interfere in such a way that only components 
with 0AM belonging to the Fibonacci sequence survive. 
As a result the pump beam entering the crystal is a su¬ 
perposition of Fibonacci-valued 0AM states, |^n)- 
After the crystal, the signal and idler states have 0AM 
values that must sum to a Fibonacci number, but are not 
necessarily Fibonacci numbers themselves. Filters in the 
signal and idler paths then remove the non-Fibonacci val¬ 
ues. The result is that at the output of the source there 
is an entangled two-photon state with Fibonacci-valued 
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0AM in the two output directions: 

V; = ^ (|F„_iU|F„_2)b + |F„_2)^|F„_i)b), (1) 

n 

where the index n runs over the indices of the al¬ 
lowed Fibonacci numbers in the pump beam: \^)pump = 
Sn 1^^)' Alice and Bob each receive half of the entan¬ 
gled pair. Each of them uses a 50/50 nonpolarizing beam 
splitter to randomly direct the photon either to one of two 
types of detection stages, referred to as L and D stages. 
L-type detection consists of an 0AM sorter [iM3 foi- 
lowed by a set of single-photon detectors. 0AM sorters 
send different I values into different outgoing directions, 
so that they register in different detectors, allowing the 
value of I to be determined. In contrast, B-type detec¬ 
tion is used to distinguish between different superposition 
states \Sn) of the form 

|5„) = /=(|F„_i) + |F„+i)). (2) 

The detection of such superposition states can be ac¬ 
complished in several ways (Sin. One complica¬ 
tion with the B-type detection is that adjacent su¬ 
perposition states, such as ^ (|F^_i)-h |F/+i)) and 

^ (l^n+l) + I F/+ 3 )) are not orthogonal, meaning that 
the two states cannot be unambiguously distinguished 
from each other. It is also necessary to keep the possi¬ 
ble key values uniformly distributed, so that Eve can’t 
obtain any advantage from knowledge of the nonuniform 
distribution. These complications add some complexity 
to the classical exchange (see [9[ for details) and alter 
the corresponding detection probabilities; however this 
extra complexity is well compensated by the increased 
key-generating capacity. Moreover, the degree of com¬ 
plication does not grow with the size of the alphabet 
used, so that the benefits outweigh the complications by 
a larger amount as the range of I values increases. 

Note that the action of the spiral does not lead to loss 
of any photons or energy. The energy is simply being 
redirected from non-Eibonacci to Eibonacci modes via 
constructive and destructive interference, causing no re¬ 
duction of photon efficiency in the apparatus. The sig¬ 
nal and idler filters after the crystal, on the other hand, 
do lower the efficiency through photon loss. The per¬ 
cent of photons retained is simply given by the frac¬ 
tion of the values in the chosen operating range that fall 
on the Eibonacci sequence; for example if the eight Ei¬ 
bonacci values between 2 and 55 are used, this about 
15%. Since no losses occur in the spiral, this is not 
enough to lower the event rate below reasonable levels. 
As in all entanglement-based QKD protocols, the princi¬ 
pal constraints on key generation rate come simply from 
the low efficiency of the down conversion process itself 
and from losses in propagation between the source and 
receivers. 

Assume that the pump spectrum is broad enough to be 
approximately flat over a sufficient span to produce sig¬ 
nal and idler 0AM values of uniform probability over 
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FIG. 2: (Color online) Apparatus for statistically distinguish¬ 
ing sets of vertically or horizontally polarized photons from 
diagonally-polarized photons. The half-wave plates rotate 
vertical and horizontal polarization states to a diagonal state, 
in order to restore indistinguishability between paths. If the 
photons are vertically or horizontally polarized, the two de¬ 
tectors C and D will register the same number of events. For 
diagonally polarized states at -|-45° from the horizontal, de¬ 
structive interference will prevent detector D from registering 
any events, with constructive interference occurring at detec¬ 
tor C. For inputs states polarized along the other diagonal 
(—45°) the roles of the detectors are reversed: constructive in¬ 
terference occurs at D and destructive interference prevents 
C from firing. 


the range Fmo to for some mo- The out¬ 

comes for L-type detection (0AM eigenstates) that will 
be used for key generation by Alice and Bob are then sim- 
ply \Fmo), \Fmo+i), ■■■, l-Fmo+iV-i). The outcomes for D- 
type detection (two-fold 0AM superposition states) used 
by Alice and Bob for key generation run from 

l^rno) = {|^mo-l) + |^mo + l)} (3) 


to 


|5'mo+Ar-l) = {\Fmo+N-2) + \Fmo+N)} • (4) 


III. DISCRIMINATING SUPERPOSITION 
STATES 

The sorting of the 0AM eigenstates |F^) in the L ba¬ 
sis is straightforward, but the sorting of the superposition 
states \Sn) in the B basis is more problematic: because 
these states are not orthogonal with the neighboring 
states two units above and below them ((S'n± 2 | 5 'n) = \) 
they can not be distinguished unambiguously from each 
other. The complications this brings are the price of in¬ 
creasing the size of the coding space. 

To deal with this, we work by analogy to superposi¬ 
tion states of polarization. Suppose Alice is sending a 
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FIG. 3: (Color online) Two copies of the previous interfer¬ 
ometer (Fig. [2]) can be combined, with the half-wave plates 
in the upper and lower copies taking the polarization vectors 
to opposite diagonals. If either of the states \V) or \H) is 
input, all four detectors have equal probability of firing. But 
if one of the two superposition states | = -^ (| 4) + | -^)) 

or I \) = (I t) “ I -^)) is input, then only the C detec¬ 
tor in one interferometer and the D detector in the other can 
fire. (The tt phase shift in one branch of the bottom interfer¬ 
ometer is not necessary, but is inserted to keep the labeling 
of the C and D detectors consistent in the upper and lower 
interferometers. ) 


stream of photons to Bob, and that each photon is either 
polarized in the vertical/horizontal basis or in the diag¬ 
onal basis. Imagine that Bob sends each photon he re¬ 
ceives through a Mach-Zehnder interferometer, as shown 
in Fig. [21 In this interferometer, a polarizing beam split¬ 
ter (PBS) transmits vertically polarized photons and re¬ 
flects horizontally polarized photons. After separating 
the two polarization components, the each component is 
rotated by a half-wave plate so that the two polariza¬ 
tion vectors point along the same diagonal. It is then 
impossible to determine which path was taken to reach 
the second (nonpolarizing) beam splitter (BS). The am¬ 
plitudes in the two arms of the interferometer are then 
recombined and sent to the two detectors labeled C and 
D. If the incoming photon was either vertically or hor¬ 
izontally polarized, then the two detectors are equally 
likely to fire. On the other hand, if the initial photon 
was diagonally polarized, | (I t) + I ~^))i then 

at the final beam splitter there will be interference be¬ 
tween the two amplitudes. In particular, there will be 
constructive interference at detector C and destructive 
interference at detector D. 

We may expand the apparatus of Fig. [2] to make 
its appearance more symmetrical and to make the anal¬ 
ogy to the 0AM case clearer below. By using a beam 
splitter to couple the interferometer of Fig. [2] to a sim¬ 
ilar one in which the half-wave plates rotate polariza¬ 
tions to the opposite diagonal (Fig. [3]), we now have a 
setup in which vertically or horizontally polarized pho- 



FIG. 4: (Golor online) A portion of a detection unit for sta¬ 
tistically detecting superpositions of 0AM states, analogous 
to the polarization version of the previous figure. The sorter 
separates different 0AM values, which are then shifted to zero 
0AM by spiral wave plates or holograms (the yellow circles). 
From the top downward, the 0AM shifters shown change the 
incoming 0AM value by A/ = —Fn- 2 , —Fh, —Fh+ 2 , —Tn+ 4 . 
Superpositions of the form \Fn) + \Fn+ 2 ) cause construc¬ 
tive interference at the Cn detectors and destructive inter¬ 
ference at the Dn detectors, while 0AM eigenstates lead to 
equal detection rates at both types. Goming out of the same 
sorter is a similar arrangement (not shown) for the states 
... |Fh_i), l-Fn+i), jFh-i-a),... 


tons are equally likely to trigger any of the four detectors, 
but the diagonally polarized superposition state | = 

(I t) + I ^)) can only trigger detectors C and Dr^. 

Similarly, superposition state | \) = -^ (| ^) — | -^)) 
can only trigger detectors and D/^. On a single 
trial, it is impossible to determine the incoming polar¬ 
ization state of an individual photon, but over many tri¬ 
als the distribution of the states received by Bob can be 
built up and compared to the distribution sent by Alice. 
This setup can therefore be used to statistically detect 
tampering with the states en route. 

A similar idea can be used for 0AM (Fig. |4]). An 0AM 
sorter replaces the polarizing beam splitter, and a pair of 
detectors Cn and Dn is used at the output ports of the 
final nonpolarizing beam splitters. If Cn fires during the 
key-generating trials, we count that as an |5'n-i) detec¬ 
tion. Due to destructive interference, Dn should not fire 
for superposition state input of the considered form, so its 
hringwill count as an \Fn) detection. Then the scheme of 
ref. 0 is used to reconcile Alice’s and Bob’s trials by clas¬ 
sical information exchange in order to arrive at an unam¬ 
biguously agreed-upon key. During the security checks, 
the distribution of counts in Cn and Dn separately are 
examined in order to detect eavesdropper-induced devi¬ 
ations from the expected probability distributions. In 
order to achieve the indistinguishability required for in¬ 
terference, the 0AM of each photon is shifted to zero 
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after the sorting (by means of a spiral phase plate, for 
example). This is analogous to the use of a diagonal po¬ 
larizer to restore indistinguishability in the polarization 
case. Note that measurements with detectors Cn and 
Dn^ respectively, are equivalent to looking for nonzero 
projections onto the states 

\Cn) = ^{\Fn) + \F^- 2 )) = i\Sn-l) (5) 

\D„) = ^(|F„)-|F„_2)). (6) 

These two sets of states are also mutually nonorthogo- 
nal; for equal n, we find {Cn\Dn) = 0, but more gen- 
erally (C„|Dm) = i { 6 m,n -2 + Sm,n+ 2 ) ■ Note that |C„) 
and \Dn) are the analogs of diagonal polarization states 
in the two-dimensional subspace spanned by |F^) and 
\Fn-2). 

The configuration of Fig. 0] is analogous to multiple 
copies of the interferometer of Fig. [2l all being fed by 
a single 0AM sorter. This generalizes Fig. [3] from two 
interferometers to many. The sorter plays the role of 
the polarizing beam splitter in Fig. [21 directing differ¬ 
ent input states into different paths through the system. 
Similarly, the diagonal polarizers and the 0AM shifters 
play the same role in the two devices, each being used 
to restore indistinguishability to the amplitudes that fol¬ 
lowed different paths, allowing these amplitudes to inter¬ 
fere when recombined at the detectors. Eigenstates are 
equally likely to trigger the C and D detectors, while su¬ 
perposition states of the form of |5'n) should never trigger 
the D detectors. As in the polarization case, the state of 
an individual trial can’t be determined, but the statistical 
distributions of detections are changed for different input 
states; the sole exception to this is that if a T) detector 
fires we know that the state must be an eigenstate. 


IV. GENERALIZATIONS 

The apparatus of Fig. |4]can be generalized in a num¬ 
ber of ways. For example, consider the setup in Fig. |5l 
Again, the 0AM sorter sends eigenstates outward in a 
series of spokes, but now a set of 50/50 beam splitters 
causes each spoke to intersect with two others. Although 
not drawn in, it is implied that there are beam splitters 
at each of the points where lines in the diagram split 
or cross. Each of the yellow circles shifts the 0AM to 
zero, so that interference can occur. Suppose the lines 
going out from the beam splitter in the portion shown 
carry angular momenta In, • • •, ^n+ 4 - Assume that 
the state entering the sorter is a uniform superposition 
of these states ^ |/n)- Then the state arriving at de¬ 

tector Dn is proportional to \ln) — \ln-i) — \ln- 2 ), while 
Cn receives i (|/n) + \ln-i) + |^n- 2 ))- If the states enter¬ 
ing the 0AM sorter are filtered to allow only a subset of 
the so-called Tribonacci numbers (obeying the three-term 
Tribonacci relation Tn = Tn-i Tn -2 + [H, 



FIG. 5: Portion of a detector setup for discriminating Tri¬ 
bonacci states. Each yellow circle represents an 0AM shifter 
that takes the 0AM in that branch to zero. Each place where 
lines split or cross is implied to have a beam splitter. The 
portion shown has 0AM values from G to /n +4 coming out of 
the 0AM sorter. 


then the setup, when fed with a source of entangled 
three-photon states, can be used to construct a three 
photon protocol based on the Tribonacci numbers, in di¬ 
rect analogy to the two-photon Eibonacci protocol. Al¬ 
though rapidly becoming less practical with increasing 
V, the extension to the N — bonacci recurrence relation, 
Tn — Tn—1 T Tn—2 T • • • Tn—N 1^ obvioUS. 

More generally still, a similar setup can be defined in 
which beam splitters cause each outgoing beam to inter¬ 
sect with any n — 1 others before reaching a detector, 
including intersections between beams that are not nec¬ 
essarily nearest neighbors. By further adding appropri¬ 
ate phase shifts and attenuation factors in the beams, we 
may use the setup to detect other n-fold linear combi¬ 
nations of 0AM states. Suppose, for example, that one 
wishes to test for linear combinations of 0AM values of 

the form ai\xm-i) + a 2 \xm- 2 ) H-h an\xm-n), where 

{xm} is some predetermined set of 0AM values. A simi¬ 
lar recursive tree can be constructed such that: (i) Each 
outgoing line intersects n beam splitters, (ii) Before the 
jth beam splitter in any line, the intensity is attenuated 
by a factor of tj = | | (where Omax is the largest of 

the coefficients) and phase shifted by 0 = arg{aj). The 
result is that the C-type detectors will be sensitive to 
superposition states of the desired form. If the superpo¬ 
sitions are non-orthogonal, then there is of course still 
ambiguity in identifying them, but over multiple trials it 
will be possible to know whether the incoming states be¬ 
long to this set of superpositions or not by their statistics, 
as in section m 

Any desired linear combination of 0AM states can be 
detected in a similar manner. As one particular exam¬ 
ple, Eig. [6] shows a portion of a setup in which the Cn 
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FIG. 6: (Color online) Portion of setup for detecting states 
of the form —\xn) + ^\xn-i): where {xn} is a predetermined 
set of 0AM values. One output from the first beam splitter 
in the nth line is attenuated in amplitude by | and the other 
output is phase shifted by tt. 


detectors test for states of the form —\xn) + 
for some known set {x^} of 0AM values. Again, if the 
input is an 0AM eigenstate with / = then the Cn 
and Dn detectors will fire with equal probability, while 
Cn+i and Dn+i will fire with one quarter of the proba¬ 
bility. In contrast, if the input state is proportional to 
— \xn)-\-\\xn-i)^ then the D detectors will not fire due to 
destructive interference. The same setup, but with the 
beam splitters mixing next-nearest neighbors instead of 
nearest neighbor beams would similarly test for states of 
the form -\xn) + ^\xn- 2 )^ 

Further note that if the additional ingredient of coin¬ 
cidence counting is added, then sufficient numbers of de¬ 
tectors and beam splitters will allow detection of various 
two-particle states in a similar manner. A two-photon 
bilinear state such as |/n+4) * (|^n) + |^n+i) + \^n-\- 2 ))^ for 
example, can be detected by connecting a coincidence 
counter between the output of detector Cn in Fig. [5] and 
a detector placed directly in the /n+4 output of the 0AM 
sorter. 

One particular example of the usefulness of this ap¬ 
proach can be seen by considering the switching between 
pairs of mutually orthogonal bases of dimension greater 
than two. In Fig. [71 detection stages for a pair of mu¬ 
tually unbiased bases of dimension four are shown. Fig. 
[7| (a) shows the setup for detection in the C basis (0AM 
eigenbasis), for the four-dimensional space spanned by 
states |F„), |F„+i), |F„+2), |F„+3). Fig. [ 7 ] (b) shows an 
arrangement that detects states in a basis that is mutu¬ 
ally unbiased with respect to jC. In this second arrange¬ 
ment, the detectors Mi, M2, M3, M4, respectively, detect 


the states 



1^1) = 

2 { + ^n+l) + ^n+ 2 ) + ^n+ 3 } 

(7) 

1^2) = 

I {|F„) + |F„+i) - |F„+ 2 ) - |F„+ 3 } 

(8) 

1^3) = 

I {|F„) - |F„+i) - |F„+ 2 ) + |F„+ 3 } 

(9) 

1^4) = 

I {|F„)-|F„+i) + |F„+ 2 )-|F„+ 3 }. 

(10) 


Random switching between the two bases can be done 
passively, by having a beam splitter randomly send each 
photon to one arrangement or the other. This is in con¬ 
trast to the setup in the original experiment of , which 
switched between two mutually unbiased bases of dimen¬ 
sion three by means of a complicated arrangement of me¬ 
chanically shifted holograms. 

V. STATE SYNTHESIS 

We note that the situation presented in the previ¬ 
ous sections contains two different sequences of linear 
combinations. There is the recursive sequence of val¬ 
ues Fn = Fn -1 Fn -2 on the space of real numbers, 
and there is the nonrecursive sequence of states \Fn) (or 
equivalently, of \Cn) and \Dn)) on the Hilbert space of 
the system. This leads to the question of whether similar 
optical arrangements to those shown above can produce 
more general linear combinations on the Hilbert space, 
and whether we can in fact construct recursive sequences 
of states in Hilbert space. In general this would mean 
a sequence of states in which each is proportional to a 
linear combination of several of the previous states in 
the sequence, according to some regular rule. We require 
proportionality rather than equality due to the need to 
normalize the states. 

With some slight alterations, the interferometric trees 
of the previous sections can be used to prepare de¬ 
sired superposition states from incoming 0AM eigen¬ 
states, instead of detecting different types of pre-existing 
superpositions. By removing the detectors and 0AM 
shifters (the yellow circles in the figures of the previ¬ 
ous sections) and by illuminating the sorter with a uni¬ 
form superposition of 0AM states, |/), we may ar¬ 
range for different superpositions to appear at the vari¬ 
ous output lines on the right. For example, the same ar¬ 
rangements described in the last section to detect states 

Cii\lm-i) F a 2 \lni- 2 ) - ^ Cin\lm-n) cau, when detectors 

and 0AM shifters are removed, prepare these same states 
at the locations previously occupied by the C detectors. 

For example, consider Fig. 0] again, but now with 
these changes. The exit ports where the C-type detec¬ 
tors previously were will now output states proportional 
to \ln) F \ln-2)- Similarly, removing the T^-type detec¬ 
tors, the corresponding output ports will produce states 
proportional to \ln) — |/n- 2 )- By building up more compli¬ 
cated trees, we can construct outputs that are arbitrary 
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FIG. 7: (Color online) (a) An 0AM sorter followed by detec¬ 
tors makes measurements in a four-dimensional space spanned 
by the four incoming basis states (the states |Fn), • • • l^hi+s) 
in the case drawn), (b) An arrangement for measuring in a 
basis that is mutually unbiased to the basis of (a). The yellow 
circles shift the 0AM to zero, so that interference can occur, 
and 50/50 beam splitters are implied at each of the points 
where lines cross. Additional beam splitters (not shown) are 
needed in some of the lines to equalize the detection proba¬ 
bilities in the four detectors (Mi,..., M 4 ). 


linear combinations of the 0AM states. This ability to 
tailor different superpositions of 0AM states may be use¬ 
ful for a number of applications in quantum communica¬ 
tion and quantum computing. 

Finally, consider the possibility of creating states that 
themselves obey a recurrence relation. As a concrete 
example, consider again the Fibonacci relation. But 
now we wish the states themeselves, not their 0AM val¬ 
ues, to obey the relation. In other words, we desire 
a set of states \xn) satisfying \xn) ^ \xn-i) + l^n- 2 ), 
where ^ denotes equality up to normalization and the 
Xn are the allowed 0AM values of the states. Suppose 
we again illuminate an 0AM sorter with a uniform su¬ 
perposition of 0AM states, |/). We then take two of 
the output lines from the sorter (say those for I = Xn 
and I = Xn-\-i) and feed them into the unit cell shown 
in Fig. [HI The input state to the unit cell is then 


INPUT 


OUTPUT 



FIG. 8: (Golor online) The basic unit cell for the recursive 
state generator (RSG) in the case of the Fibonacci recurrence 
relation. Given a superposition of two input state (for exam¬ 
ple \xn) = \Fn) and |xn+i) = |Tn+i)), the cell produces a 
state proportional to \xn-\-2) = (l^n) + |^Cn+i).)- The at¬ 

tenuation of the amplitude in the bottom line by is there 
to ensure that the amplitudes of the three states are equal. 
By adding additional attenuations and phases shifts in the 
various branches, units cells for any other complex, linear re¬ 
cursion relation can be constructed. 



Attenuation 


FIG. 9: (Golor online) Recursive state generator (RSG). Mul¬ 
tiple copies of the unit cell of Fig. [8] (represented by the 
blue boxes) can be concatenated in a fractal-like manner to 
form a set of output states that obey the Fibonacci relation 
\xn) = -^\xn-i)-\-\xn- 2 )- Each additional ccll added has half 
the output intensity of the previous cell. The attenuations 
(beam splitters), represented by the red circles, are added to 
make all of the hnal output intensities equal. Adding addi¬ 
tional attenuations and phases shifts inside the unit cell, other 
recurrence relations can be realized in Hilbert space. 


{\xn)i F 1 ^ 71 + 1 ) 2 ), where the labels 1 and 2 refer to 
the two input ports. Then in addition to the original 
superposition leaving from the top two output ports on 
the right of the cell, there is an equal amplitude for 
the superposition state lx^+ 2 ) = (l^n) F |x^+i)) to 

exit the lower port. In other words, the input state 
(l^n)i + 1 ^ 77 + 1 ) 2 ) is transformed into the output state 

{\xn)i + \xn+i )2 + \xn+ 2 )z)- Now suppose we nest 
multiple copies of these unit cells together as in Fig. El 
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FIG. 10: (Color online) The apparatus of Fig. [H (with the 
detectors removed) is abbreviated as a single unit, represented 
by the green box (left side of hgure). Multiple copies of this 
unit can be repeated, with the output of each used as input 
to the next (right side). The result is that the apparatus 
alternately switches between the L and D bases. 


If the attenuation factors are appropriately adjusted, the 
states at the output ports (assuming p unit cells are used) 
will be |x„)i, |a:;„+i) 2 , |a:„+ 2)3 •. • |a:;„+p+i)p+ 2 , each ap- 
pearing with equal amplitude. We refer to this arrange¬ 
ment as a recursive state generator (RSG). (Recursively 
constructed quantum states have been considered before 
from a different perspective.) 

As the number of unit cells nested in the RSG in¬ 
creases, the output intensity decreases exponentially: 
each additional unit drops the output per output port 
by a factor of two. But there is nothing in the setup 
that is especially sensitive to high intensities, so by using 
large input intensities it can be arranged to have a suffi¬ 
ciently large number of cells to produce complicated sets 
of states while maintaining useable output levels. 

By placing the RSG at the output of an 0AM sorter 
and allowing switching between different output lines of 
the sorter, different sequences obeying the same recur¬ 
rence relation are produced, for example switching be¬ 
tween the Lucas and Fibonacci sequences. More gener¬ 
ally, by nesting multiple recursive trees inside each other, 
this gives access to finite approximations of fractal states 
in Hilbert space. The construction of such recursive 
chains of states allows the possibility of carrying out in¬ 
formation processing or other tasks on these states in a 
relatively simple manner. 

To give one application of the approach presented here, 
consider taking multiple copies of the apparatus from Fig. 
[T] and feeding the output of each copy into the input of 
another copy in a fractal manner as shown in Fig. [TOl If 
Alice and Bob both feed their half of a down conversion 
pair into such an apparatus, the result is an implementa¬ 
tion of an entangled two-photon quantum walk in 0AM 
space: the 0AM values of each photon can walk up and 



FIG. 11: (Golor online) The apparatus of Fig. [TO] implements 
a random walk. By alternating L and D measurements, each 
of the two entangled photons can be made to walk randomly 
up and down one of the two chains of states shown. The left- 
hand chain is constructed from states \Fn) with odd index n, 
the right-hand chain is built from states with even n. The 
\Cn) and \Dn) states are those that are detected by the Cn 
and Dn detectors in Fig. [H 


down the two chains shown in Fig. [TTJ The 0AM value 
can be thought of as the walk variable, with the interfer¬ 
ometers that choose between the C and D states acting 
as the coin operators. (A quantum walk in 0AM space 
has recently been implemented by other means in |25|.) 


VI. CONCLUSIONS 


In this paper, we have introduced an interferomet¬ 
ric method for statistically distinguishing 0AM super¬ 
position states both from eigenstates and from other 
nonorthogonal superpositions. In addition, the same es¬ 
sential method can be used to prepare superpositions of 
optical 0AM states. This ability to tailor different su¬ 
perpositions of 0AM states may be useful for a number 
of applications in metrology and quantum computing; in 
particular, it allows complex sets of states to prepared in 
a simple manner, allowing the possibility of carrying out 
sophisticated information processing algorithms directly 
on the Hilbert space of the system. Further, by replac¬ 
ing the 0AM sorters with devices that sort according to 
other degrees of freedom, similar arrangements can be 
used for superpositions of eigenstates of other operators. 
The systems described can all be placed on integrated 
optical chips with current technology or with that which 
will be available soon, allowing them to be used to gener- 
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